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NOTATION 


A_  Coefficients  of  simultaneous  equations  defined  in  Equation  (36) 

B  Expression  defined  by  Equation  (13) 

Frictional  drag  coefficient 

2 

C  Drag  coefficient;  drag  divided  by  1/2  p  V 

s 

CDc  Cavity  drag  coefficients 

C  Power  coefficient  =  wQ/ (Tr/2pv\^) 

P  a 

C^,  Thrust  coefficient  =  T/(ir/2pV^R^) 

c  Chord  length  divided  by  propeller  radius 

c  Camber  correction  factor 

o 

D  Propeller  diameter 

EAR  Expanded  area  ratio 

H  Expression  defined  by  Equation  (14) 

Quantities  defined  by  Equation  (20-1) 

J  Advance  coefficient:  V^/(nD) 

2  5 

Kq  Torque  coefficient:  Q/(pn  D  ) 

,  2  4 

Thrust  coefficient:  T/ (pn  D  ) 

i  Cavity  length  divided  by  the  chord  length 

n  Revolutions  per  second,  rps 

P  Pressure 

P/D  Pitch-diameter  ratio 

2 

p  •  Nondimensional  perturbation  pressure  (P-P  )/l/2pV  ) 
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ABSTRACT 


A  lifting-surface  theory  and  numerical  procedure  for 
designing  supercavitating  propellers  are  presented. 

Both  a  subcavitating  and  a  supercavitating  propeller 
are  represented  by  vortex  and  source  distributions.  Unlike 
the  subcavitating  propeller,  however,  source  strengths  for 
a  supercavitating  propeller  are  related  to  cavity  thick¬ 
ness,  which  is  not  known  without  examination,  and  have  to 
be  obtained  by  solving  related  integral  equations. 

Numerical  solution  of  the  integral  equations  is  obtained  as 
a  correctional  function  for  a  stripwise  supercavitating 
cascade  theory  which,  with  lifting-line  theory,  is  used  for 
preliminary  design  of  supercavitating  propellers.  The  in¬ 
duced  axial,  radial,  and  tangential  velocities  are  obtained 
on  a  blade  reference  surface  that  allows  arbitrary  skew, 
rake,  and  radial  pitch  variations.  The  blade  shape  is 
obtained  -as  a  correction  to  the  shape  obtained  from  strip- 
wise  supercavitating  cascade  theory.  Thrust  and  torque 
coefficients  are  obtained  from  pressures  on  the  blade 
surface. 

The  method  is  applied  in  designing  several  super¬ 
cavitating  propellers  that  have  design  conditions  the  same 
as  those  of  existing  supercavitating  propellers. 

Numerical  computations  were  also  used  to  design  two 
additional  supercavitating  propellers  which  were  built  and 
tested.  The  design  predictions  are  compared  to  the  experi¬ 
mental  data,  both  for  blade  cavity  height  and  performance, 
characteristics,  and  good  correlation  is  obtained. 

ADMINISTRATIVE  INFORMATION 

The.  Naval  Material  Command  (NAVMAT  08T)  funded  this  investigation  in  support 
of  an  ongoing  Ship  Performance  and  Hydromechanics  Exploratory  Development  Program 
(Program  Element  62543N,  Task  Area  ZF43421001)  assigned  to  the  David  W.  Taylor  Naval 
Ship  Research  and  Development  Center  (the  Center) .  This  work  has  been  completed 
under  the  High-Speed  Propulsor  Task  portion  of  this  program,  (Center  Work  Unit 
1500-103),  and  partly  under  the  Naval  Sea  Systems  Command's  General  Hydrodynamic 
Research  Program  (Center  Work  Unit  1542-817). 

INTRODUCTION 

Ever  since  supercavitating  propellers  were  first  investigated  systematical- 

12*  3  4 

ly,  *  steady  progress  has  been  made  in  understanding  associated  problems.  ' 

Designers  of  supercavitating  propellers  have  been  urged  to  consider  the  interference 


*A  complete  listing  of  references  is  given  on  page  83. 
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effect  of  neighboring  cavities  and  the  effect  of  flow  retardation.  ’  To  this  end, 
three-dimensional  cavity-flow  theory  similar  to  the  lifting-surface  theory  for  sub- 
eavitating  propellers'’’  has  been  formulated. ^  However,  to  date  none  of  the 
theories  has  been  applied  to  the  actual  design  of  a  supercavitating  propeller  owing 
to  the  numerical  complexity  of  the  problem. 

As  with  subcavitating  propellers,"’  a  supercavitating  propeller  may  be  designed 

in  two  steps:  preliminary  or  lifting-line  design  and  final  or  lifting-surface 

design.  The  former  gives  an  approximate  solution  and  supplies  basic  data  to  the 

latter.  The  preliminary  design  supplies  information  for  the, final  design  concerning 

load  distribution  on  the  blade  and  preliminary  pitch  distribution,  which  forms  the 

basic  singularity  surface  for  the  final  design.  The  section  shape  of  the  blade 

cavity  and  the  cavity  length  of  a  supercavitating  cascade  model  are  also  computed  in 

the  preliminary  design.  The  information  is  used  in  the  final  design  as  the  first 

approximation.  Although  this  report  deals  mainly  with  the  final  design,  we  also 

11 

consider  aspects  of  the  preliminary  design  to  help  the  reader  understand  the  final 
design  more  fully. 

The  propeller  diameter,  the  blade  number  and  contour,  the  hub  diameter  for  the 
given  propeller-thrust  power  coefficient  with  the  given  ship  speed,  propeller  revo¬ 
lutions  per  minute  (rpm) ,  and  the  wake  fraction  are  considered  to  have  been  deter¬ 
mined  before  applying  the  present  preliminary  and  final  programs.  To  ensure  that 
the  designed  propeller  is  fully  cavitating  and  that  the  blades  are  thick  enough,  the 
leading-edge  cavity  thickness  is  specified.  If  the  specified  cavity  thickness  near 
the  leading  edge  does  not  accommodate  a  stable  cavity  of  at  least  1.5  chord  lengths 
because  the  cavitation  number  is  too  large,  then  the  leading-edge  cavity  thickness 
is  increased  internally  to  make  a  long  enough  cavity  for  the  preliminary  program. 

The  purpose  of  the  present  program  is  to  design  a  fully  cavitating  propeller 
that  is  efficient,  has  no  face  cavitation,  is  structurally  strong,  and  meets  the 
design  requirements. 

As  for  a  subcavitating  propeller,  the  supercavitating  propeller  is  represented 

g 

by  vortex  and  source  distributions.  The  source  strengths  are  related  to  the  cavity 
thickness,  are  not  known  a  priori  and  are  to  be  obtained  by  solving  related  singular 
integral  equations.  When  the  source  strengths  are  known,  the  supercavitating  pro¬ 
peller  problem  is  similar  to  the  problem  for  a  subcavitating  propeller  with  thick 
rr.d  wide  blades.  Thus,  many  of  the  computation  techniques  apply  to  both 


.2 


subcavitating  5  and  supercavitating  propellers.  Indeed,  this  present  study  makes 
use  of  many  parts  of  the  lifting-surface  design  programs  reported  by  Kerwin  for 
a  subcavitating  propeller.  The  coordinate  systems  of  the  blade  are  the  same  for 
both  propeller  types. 

Since  the  singular  integral  equation  is  a  Fredholm  equation  of  the  first  kind, 
the  method  for  its  solution  must  be  chosen  with  extreme  care.  Additionally,  the 
cavity  shape  is  not  known  without  examination.  Thus,  iteration  and/or  some  special 
cavity  model  has  to  be  applied  knowing  that  all  the  inviscid  cavity  models  are  not 
exact;  rather,  they  are  approximate  representations.  The  source  distribution  is 
obtained  from  a  stripwise,  two-dimensional,  supercavitating-cascade  representation 
developed  in  the  preliminary  design.  For  the  procedure  used  at  present,  the  source 
distribution  is  multiplied  by  a  double  polynomial  having  unknown  coefficients  in 
terms  of  the  radial  and  chordwise  coordinates.  Since  the  two-dimensional  cavity 
model  used  in  preliminary  design  is  a  linear,  double  spiral  vortex  model,  the  cavity 
streamline  is  closed  at  infinity  instead  of  at  the  cavity  end.  Thus  the  unknown 
source  strengths  are  distributed  in  a  prefixed  plane  that  contains  the  blade  surface 
and  extends  to  the  wake  about  1-1/2  chord  lengths.  At  the  end  of  the  plane  is 
attached  another  source  line  having  an  unknown  polynomial  strength  that  is  to  be 
solved  along  with  the  double  polynomial  coefficients  by  the  least  squares  method. 

The  solution  is  a  function  of  propeller  geometry,  given  thrust  or  power,  blade  load 
distribution,  advance  coefficient,  and  cavitation  number. 

The  induced  axial,  radial,  and  tangential  velocities — thus  pitch  and  camber 
distribution— are  obtained  on  a  blade-reference  surface  that  allows  arbitrary  skew, 
rake,  and  radial  pitch  variation.  The  blade  cavity  shape  is  obtained  as  a  correction 
to  the  blade  cavity  derived  from  supercavitating  cascade  theory.  The  thrust  and 
torque  coefficients  are  obtained  from  the  pressure  on  the  blade  surface,  which  is 
converted  from  the  lift  distribution  through  the  Kutta-Joukowsky  theorem. 

If  the  hub  boundary  condition  is  not  considered  in  solving  for  the  cavity 
source  distribution,  the  radial  velocity  caused  by  the  cavity  source  is  not  stable. 
When  the  hub  boundary  condition  is  considered,  the  numerical  results  are  shown  to 
have  reasonable  convergence  in  many  numerical  experiments.  Two  new  supercavitating 
propellers  were  designed  using  numerical  computations  derived  during  the  present 
program,  and  models  were  manufactured.  Model  test  results  correlate  reasonably  well 
with  the  computed  cavity  shapes  and  powering  performance. 
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GEOMETRY  OF  THE  BLADE 


The  geometry  of  the  blade  of  a  supercavitating  propeller  is  almost  the  same  as 
that  of  a  subcavitating  propeller.  The  blade  is  represented  in  the  flow  by  singu¬ 
larity  distributions  of  vortices  and  sources  on  the  blade-reference  surface  in  the 
linear  version.  Thus,  the  blade-reference  surface  has  to  be  known  approximately, 
although  the  blade  shape  has  to  be  obtained  as  a  solution.  As  in  the  theory  for  sub¬ 
cavitating  propellers,  the  reference  surface  will  be  close  to  the  helical  surface  of 
the  hydrodynamic  pitch  angle  3^  obtained  in  the  preliminary  design.  Although  the 
reference  surface  could  be  determined  more  accurately  by  an  iterative  procedure,  in 
the  present  program  3^  Is  used  once  without  iteration,  which  is  considered  to  be  a 
reasonable  approximation.  The  pitch  angle  of  the  reference  surface  of  the  wake  can 
be  set  differently  from  that  of  the  blade. 

A  cylindrical  coordinate  system  (x,r,0),  seen  in  Figure  1,  is  defined;  the  x- 
axis  is  coincident  with  the  axis  of  rotation  of  the  propeller  and  is  positive  when 
facing  downstream.  Thus,  the  reference  surface  can  be  represented  by 

x  =  A0  +  ?( r),  rR  <  r  <  1,  0L  <  6  <  00  (1) 

A  =  i'  tan  3,  (2) 

1 

where  £(r)  depends  on  rake  and  skew,  and  A  is  in  general  close  to  a  constant.  If  a 
Cartesian  coordinate  is  defined  as  in  Figure  1,  the  parametric  representation  of  the 
reference  surface  x(r,0)  =  xi  +  yj  +  zic  is 

x(r,0)  =  A0  +  £(r) 

y(r,0)  =  r  cos  0  (3) 

z(r,0)  =  r  sin  0 


Then  the  area  element  is 
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where 


„  ~  -  2  2  q 

E  =  Xg  •  Xg  =  r  sec  0^ 

G  =  x  •  x  =  (6X  +£  )2  +  1 
TC  x  v  r 

F  =  x  •  xa  =  X(0X  +C  ) 
r  0  r  r 


Thus 


1/2 

dS  =  [X2+r2+r2(0Xr+Cr)2]  drd0 

The  unit  vector  normal  to  the  reference  surface  is 

Xg  x  Xg  x  xr  (-ri+r(Xr0+£r)j+Xk 
lx0Xxri  /eG-F2 

If  X  and  £  are  independent  of  r,  the  normal  vector  in  Equation  (7)  is  on  the  cyl¬ 
inder,  r  =  constant.  That  is,  the  radial  component  of  the  normal  vector  is  zero 
Two  convenient  unit  vectors  tangential  to  the  reference  surface  are 


(9) 


x 

r 


/G 


The  curvilinear  distance  is  measured  along  the  helix  on  the  reference  surface 
from  the  0=0  plane  to  the  leading  and  the  trailing  edges,  denoted  by  S.  and  S  , 

X/  L 

respectively.  The  chord  length  of  the  propeller  blade  at  any  radius  is 


c(r)  =  St(r)  -  S^(r) 


(10) 


VELOCITIES  DUE  TO  SINGULARITIES 

Both  a  source  distribution  with  strength  m(r,0)  per  unit  area  and  a  vortex 
distribution  with  strength  G(r,0)  per  unit  area  are  distributed  on  each  blade- 
reference  surface.  Then,  the  velocity  um  caused  by  the  source  distribution  is 
derived  from  potential  theory 


um(x,r,0)  = 


(ID 


— G 

The  velocity  u 
Savart  law 


caused  by  the  vortex  distribution  is  derived  from  the  Biot- 


uG(x,r,0) 


where 


B  =  j(x-£)2+r^+p2-2rp  cos(<|)+6^-0) 


1/2 


«  2k(k-l) 

K  Z 


(12) 


(13) 


(14) 


(15) 


I HMM  lar.wrvrtrrwr  |  i:rr*nyw^ 


Vi 

'9  * 


B  =  {x-£(p)}i  +  {r-p  cos  (4>-HS,  — 0 ) }j  -  (p  sin  (<j)-6, -0)}k 


where  Z  is  the  number  of  blades 


d£  =  {tan  6i+  cos  ($+6^-0) j+  sin  (<J>+6k~0)k}dp 


which  is  a  line  element  along  the  intersections  of  the  0  =  <p  plane  and  the  reference 
surface.  6  is  the  angle  between  d£  and  the  x  =  constant  plane  or  tan  6  =  9x/9r  = 

+  £  (r).  The  axial,  tangential,  and  radial  components,  u  ,  u„,  and  u  are  written  as 

1“  Si  t/  X 

follows 


1  00 


u  (x,r ,0) 

a 


‘hi!  »(p.«273 (x-5> h 


rH  6L 


1  00 


Uq (x,r,0)  ~  ~  ^  j  J  ^(P^)  ^3  sin  (<!)+5k-0)  H  d(f>dp 


rH  0L 


ur(x,r,0) 


X  w 

=  ^TT  I  /  m(p,4>)  ^3  {r-p  cos  (<j>+6k-0)}  H  d<J)dp 


rH  6L 


1  0T 


u  (x,r,0)  = 

a 


~hj  i  G(Pi'»  fr  sin  (!!+(ik"f))i  H  d'I'dp 


rH  °L 


1  00 


7—  f  (  G  (p,$)  ^ ^  (p2-pr  cos  (4>+6, -0 ) }  H  d<j>dp  (18) 
47T  ^  K  S  “  B3  k  (cont.) 


rH  0L 


7 


1  6T 


Uq (x,r,0) 


=  Tv  j  J  G(p>$)^^3  [{r  tan  Mx-S 


rH  0L 


-*  _ %>  _»  *  ■«  „.«e 


x  cos  (4>-hS ^-0 ) } -p  tan  <S  cos  ($+6^.-6)]  H  d<f>dp 


1  «  Z 

"  I  I  Ip  tan  cos  W^k”8^ 

rH  «L  k=lB 


-  (x-o  sin  (<J)+6^-0)-r  tan  3^]  H  d<j>dp 


1  0T 


ur(x,r,0) 


■  Ij-  J  |  <XM)  £^5  {(x-5)  sin  «+6k-e 


rH  6L 


+  p  tan  6  sin  (<}>+<$  ^-9))  H  d<j>dp 


1  00 


h  f  f  °s<-l>A)J£b tp  tan  ei sin  W+Ve) 


rH  \ 


+  (x-£)  cos  ($+<5^-0)}  H  d<j)dp 


where  ru  =  hub  radius 
ri 

Gg  =  -  dG(p,<j>)/dp  =  trailing  vortex  strength  in  @L  <  <f>  <  00 

subscripts  a,  0,  and  r  indicate  axial,  tangential,  and  radial  components, 
respectively 

The  induced  velocity  component  due  to  the  singularity  distributions,  normal  and 
tangential  to  the  blade-reference  surface,  can  then  be  written  in  a  linear  approxi¬ 
mation  from  Equations  (5)  and  (7) 
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(19) 


u  = 


n  H(r,9) 


rvX(r)vr 


(9 


dX  * 

dr  dr  / 


u 

r 


u 

r 


1 

Hjfr.e) 


(XVrue) 


where 


=  r  sec  gi 


(20) 


(20-1) 


CAVITY  BOUNDARY  CONDITION 

The  cavitation  number  defined  with  respect  to  ship  speed  Vg  is  written 


°,= 


(21) 


where  is  the  pressure  infinitely  far  upstream  and  Pc  is  the  pressure  on  the 
cavity.  From  the  Kutta-Joukowsky  theorem  the  vortex  distribution  can  be  written  in 
a  nondimens ional  form  (Appendix  A)  as 


G_ 

V 

s 


1  ,  .  ,  jV_ 

2  V  Vn 


where 


P  = 


2  2, 

V  =  (V  +r  to  ) 
a 


1/2 


(22) 


(23) 
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P  is  the  pressure  on  the  foil,  and  is  the  local  velocity.  In  designing  a  foil, 
we  cr..i  use  the  local  velocity  obtained  during  preliminary  design 


V*  =  {(Va+Ua}  +(raHV  } 


1/2 


(24) 


where  V  =  V  (1-w) 
a  s 

w  =  wake  fraction 
0)  =  angular  speed 

When  the  supercavitating  propeller,  advancing  with  velocity  V  ,  is  represented 

SL 

by  source  and  vortex  distributions,  the  perturbation  velocity  component  parallel  to 
the  blade-reference  surface  can  be  written  by  a  linear  approximation  of  the  Bernoulli 
equation  of  flow,  referring  to  a  moving  frame  of  reference  (Appendix  A),  as 


m  G 

UT  UT  £  Vs 

r  + v~  =  ■  2  v~ 

s  s 


on  the  pressure  side  of  the  blade,  and 


m  G 

V  +  V 
s  s 


als 
2  V 


on  the  cavity,  since 


G  G  _ 
UT+  "  UT-  =  G 


G  .  G  _  0 
UT+  +  UT-  "  2  Uf 


(25) 


(26) 


(27) 


where  u^  is  the  tangential  velocity  due  to  the  vortex  distribution  at  all  points 
except  the  point  concerned;  subscripts  "plus"  and  "minus"  indicate  the  values  on  the 


upper  and  lower  sides  of  the  mean  camber  of  the  blade,  respectively;  from  Equations 
(25)  through  (27)  the  tangential  velocity  due  to  the  source  distribution  alone  can 
be  written  as 


and 


V 

s 


JL_  V  a  vj. 

2  V  V  2  V 


V 

s 


on  the  blade  plane  (28) 


V  2  V 
s 


V 

s 


on  the  cavity  plane  (29) 
behind  the  blade 


When  G  and  0  are  known,  these  two  equations  will  form  a  system  of  Fredholm  integral 
equations  of  the  first  kind  for  the  source  distribution. 

When  the  propeller  blade  shape  is  given  instead  of  the  vortex  distribution,  the 
problem  becomes  one  of  prediction  rather  than  design,  and  the  boundary  condition  on 
the  pressure  side  of  the  blade  is 


u 

n 


V 

s 


dF 

dx 


(30) 


where  F  is  the  blade  ordinate  with  respect  to  the  reference  surface  and  x  =  p<j>  sec 

In  this  case  the  blade  reference  surface  must  be  in  close  agreement  with  the  mean 
camber  and  the  wake  surfaces.  Now  Equations  (28)  through  (30)  become  the  integral 
equations  for  both  sources  and  vortices.  Thus,  the  prediction  problem  becomes  more 
complicated.  In  this  report  only  the  design  problem  is  considered.  However,  the 
prediction  problem  can  be  solved  similarly. 

SOLUTION  FOR  SOURCE 

For  the  lifting-surface  design  of  supercavitating  propellers,  the  load  distri¬ 
bution  p  +  0  on  the  blade  is  supplied  from  the  preliminary  design  computations.  As 
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for  subcavitating  propellers,  the  pitch  angle  fk  of  the  blade-reference  surface  is 
taken  as  the  hydrodynamic  advance  angle  given  by  the  preliminary  design  calculations. 
Then,  were  the  source  distributions  known,  all  the  induced  velocities  could  be 
computed,  and  the  blade-cavity  shape  and  the  final  pitch  could  be  obtained  in  a 
manner  similar  to  that  used  for  a  subcavitating  propeller. 

Two  problems,  strength  and  the  location,  are  related  to  the  source  distribu¬ 
tions.  The  location  may  be  considered  to  be  on  the  blade  reference  surface.  But, 
how  far  downstream  should  the  sources  be  distributed?  A  two-dimensional  supercavi- 
tating  cascade  section  could  be  considered  to  contribute  to  the  source  distribution. 
When  the  Riabouchinsky  cavity  model  is  used,  the  cavity  source  will  be  confined 

inside  a  finite  cavity  domain.  Thus  the  cavity  domain  supplied  by  two-dimensional 
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theory  could  be  used  in  an  attempt  to  solve  the  source  strength  distribution. 

If  the  linear  double  spiral  vortex  model  is  used,  the  source  distribution  does  not 
terminate  at  the  cavity  end,  but  the  wake  source  continues  from  the  end  of  the 
cavity.  Although  there  exists  a  logarithmic  singularity  in  the  normal  velocity  at 
the  cavity  end,  the  streamline  represented  by  the  integration  of  the  normal  velocity 
continue’;  smoothly  at  the  cavity  end.  Any  cavity  mode  is  known  to  give  reasonable 
predictions  of  lift  and  drag.  Therefore,  it  may  be  better  to  free  ourselves  from  the 
concept  of  a  finite  cavity  platform  and  to  consider  the  cavity  wake-source  distribu¬ 
tion  as  beginning  behind  the  leading  edge.  The  major  question  is  where  to  truncate 
and  how  to  reduce  the  truncation  error. 

Even  if  the  problem  of  the  source  domain  were  solved,  the  method  of  solution 
would  be  far  from  definite  because  there  is  no  established  method  to  solve  the 
Fredholm  singular  integral  equation  of  the  first  kind;  there  are  suggestions,  how¬ 
ever,  that  a  series  of  eigenfunctions  could  be  used.^  Thus  solutions  to  the  three- 
dimensional  problem  of  a  supercavitating  foil  are  complicated  and  time  consuming. 

To  arrive  at  a  reasonable  method  to  solve  the  present  problem,  we  considered  a 
combination  of  three-dimensional  corrections  to  a  two-dimensional  solution  and  the 
series  of  eigenfunctions.  That  is,  the  two-dimensional  solution  for  the  supercavi¬ 
tating  cascade  is  miltiplied  by  a  double  polynomial  function  with  unknown  coeffi¬ 
cients 


J 


f  (P  ,<f>)  = 


■z 


i=l 


J  J 


i=l  j=2 


i-1  j-1 

p  xJ 


in  the  blade  plane 


(31) 
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with  x  =  H.j4>,  which  is  zero  at  the  trailing  edge*  and 

J  J  J 

f  (p,<j>)  =  an  P*  ^  ^b^  p*  1  x*  ■*■  in  the  cavity  plane  (32) 

i=l  i=l  j=2 

The  two-dimensional  solution  is  supplied  from  preliminary  design  of  each 
section.  The  cavity  plane  has  to  be  finite  numerically.  A  separate  source  distri¬ 
bution  along  the  truncation  line  of  the  cavity  is  considered  with  unknown  coeffi¬ 
cients  to  compensate  for  truncation  of  the  cavity: 

J 

f2(p)  =  HA<j)^Tai  p1'1  (33) 

i=l 


In  addition,  the  cavity  thickness  near  the  leading  edge  is  specified  to  be  equal  to 
that  computed  from  cascade  theory  during  preliminary  design.  This  thickness  is 
chosen  because  many  sets  of  cavity  thicknesses  give  the  same  load  distribution  and 
the  same  cavitation  number  for  the  linearized,  two-dimensional  cavity  problem.  The 
leading  edge  conditions  can  be  interpreted  as  follows:  the  source  distribution  at 
the  leading  edge  is  exactly  the  same  as  that  of  a  tv?o-dimensional  supercavitating 
cascade,  that  is 


f  (p » ^ ) at  tjie  ieacjing  edge  ^ 

With  the  previously  described  representations  of  the  source  distribution,  the 
integration  of  each  term  in  the  integral  equation  is  performed  numerically.  If  the 
propeller  blades  have  rake  or  skew,  £  is  a  function  of  p.  This  is  treated  by  a 
linear  approximation  in  each  Ap  interval 

5  =  c0  +  p  q  (34-d 


where 


13 


Ci  =  (^0+C)._.  =  tan  6 

1  3p  p=P1 


h  =  (X0+op=Pl  _  P1  tac  6 


If  X  =  p  tan  3i  is  assumed  to  be  a  constant  in  interval  Ap,  all  the  integration  with 
respect  to  p  can  be  performed  analytically  in  the  interval  both  on  the  blade  and  in 
the  wake. 

That  is,  by  substituting  Equation.  (34-1)  in  B 


J  B3 


dp  =  P(n+1) 


where 


B  =  (ap  +2bp+c) 


a  =  1  +  tan  6 


b  =  -  {r  cos  (<{>+<5^-0)  +  (x-CQ)  tan  6} 


c  =  r2  +  (x-£  )2 
o 


_1 _  ap+b 

,  2  B 


p(2)  =  — -1-  Mt£ 

K  J  .2  B 

ac-b 


P  (3)  =  -7 


ap+b  +  g\  +  (2b  -ac)  p+bc 
/a  )  a (ac-b2)  B 
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P(n+1)  =  A 


+  B0  P<n>  +  E  P(n”1) 


where 


A  = 


(n-2)  a 


(2n-3)  b 
(n-2)  a 


E  =  -  iSliLS. 

(n-2)  a 


n  >  3 


The  integration  with  respect  to  <f>  is  performed  using  the  trapezoidal  rule.  Then  the 

integral  equation  can  be  represented  as  a  system  of  linear,  simultaneous  equations 

for  the  unknown  coefficients  a..,  b.  ,,  and  a, 

ij  ij  i 


1  J  P  Q-  K 

x^7aiiAl,m(i-n  +  J^7bpqA£m(pq)  +  \Alm(k) 

i-1  j=l  P=1  q=l  k=l 


V  V 
_G  J.  ,  OJ 

~  2V  V  2  V 
s 


u 


m 


u„ 


_f  _  _T 
V  “  V 


s  s 

on  the  blade  plane 


(36) 


aZs 

2  V 


m 

Uf  _  UT 


s  s 
on  the  cavity  plane 


where  the  A^  are  the  integrations  of  the  double  integrals  of  source  distribution  in 
Equations  (18),  (20),  (27)— (29)  both  on  the  blade  and  the  cavity  planes  at  the  &th 
collocation  point,  and  m(ij)  is  the  index  corresponding  to  the  coefficient  a^ . 

By  the  least-squares  method  of  solving  simultaneous  equations,  a  square  matrix 
(Bpq)  is  made  out  of  the  generally  nonsquare,  matrix  (A^)  where 


pq 


Li 

£ 


An  An 
£p  &q 


(37) 
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and  the  right  side  of  the  matrix  is 


L 

Bp  =  £  AJtp 
% 

where  represents  the  right  side  of  simultaneous  Equations  (36).  The  solution  of 
the  simult.  neous  equation 

./^?ii  Bpq(ij)  Bpq(£m)  Bpq(k)  3p 

p  -  1,  2,.... 

is  the  desired  least  squares  solution  for  the  cavity  source. 

IMAGES  FOR  THE  HUB 

To  satisfy  the  hub  boundary  condition 

u  =  0  on  r  =  r„ 
r  H 

the  two  image  systems  for  vortex  and  source  distributions  were  considered  separately. 

The  hub  images  for  the  vortex  distribution  are  the  same  as  those  Kerwin8  used. 

That  is,  the  same  vortex  strength  as  that  distributed  on  the  blade  (x,r..,0)  is  used 

2 

at  the  image  point  (x1#ro,0)  where  r,  •  r0  =  ru.  Then  the  trailing  vortices  tend  to 
cancel  the  radial  velocity  on  the  hub.  This  is  an  approximation  to  a  two-dimensional 
vortex  outside  a  circular  cylinder. 

The  hub  images  for  the  source  distribution  need  special  consideration  because 

the  cavity  sources  for  a  supercavitating  propeller  are  much  stronger  than  those  of 

subcavitating  propellers.  The  first  approximation  of  the  image  source  may  be  taken 

from  an  approximation  of  the  hub  by  a  sphere.  The  image  system  for  a  point  source 

2 

of  a  unit  strength  at  r  =  r^  of  a  sphere  is  a  point  source  at  =  r^/r^  the 

strength  r^/r^  anc*  a  line  sink  stretching  from  r  =  to  the  center  of  the  sphere 
with  the  strength  l/r...  In  this  way  the  total  absolute  strength  of  the  line  sink  is 
equal  to  the  point  source  at  r  =  x^.  Thus,  on  the  sphere,  r  =  rH>  the  normal 
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velocity  induced  by  the  source  outside  of  the  sphere  is  canceled  by  the  image  source 
system;  however,  outside  the  sphere,  the  velocity  field  due  to  the  image  system 
decays  as  flow  caused  by  a  doublet  inside  the  sphere. 

The  image  source  distribution  is  multiplied  by  the  same  double  polynomial  as 
Equations  (31)  and  (32),  and  the  hub  boundary  conditions  are  considered  together  with 
the  cavity  boundary  conditions,  Equations  (28)  and  (29),  to  solve  the  simultaneous 
Equations  (33),  (34),  and  (36)  by  the  least-squares  method.  The  solution  is  checked 
to  see  if  it  actually  satisfies  the  hub  boundary  condition  in  addition  to  the  cavity 
boundary  conditions. 


BLADE  SECTION  SHAPE 

When  the  unknown  coefficients  are  obtained  from  Equation  (38),  the  source 
distribution  is  computed  from  Equations  (31)  and  (33).  Thus,  all  the  induced  veloc¬ 
ities  can  be  calculated.  It  is  possible  to  obtain  the  section  shape  of  the  propeller 
blade  by  integrating  the  velocity  field  along  the  blade  reference  surface.  However, 
this  may  require  a  large  amount  of  velocity  information.  Thus,  the  method  used  by 

fL 

Kerwin  is  followed  to  correct  the  foil  shape  derived  from  supercavitating  cascade 
theory  in  the  preliminary  design  process.  At  a  field  point  of  each  section  the 
normal  velocity  component  V^n  is  obtained  on  the  pressure  surface  of  the  blade  from 
the  lifting-surface  theory,  and  V^n  is  matched  with  the  corresponding  normal  velocity 
obtained  from  supercavitating  ^ascade  theory 

m 

5^ a.  x*  +  c  v0  =  v„  (39) 

^ i.  o  2n  3n 

i=o 


where  a.  and  c  =  unknown  coefficients 
1  o 

v„n  =  normal  velocity  component  on  a  foil  of  a 
supercavitating  cascade 

x  =  distance  from  the  leading  edge 

The  unknown  coefficients  are  determined  by  the  least-squares  method.  That  is,  the 

coefficient  c  and  the  terms  a.  x1  are  the  correction  to  the  normal  velocity  compo- 
o  x 

nent  of  the  two-dimensional  supercavitating  cascade  caused  by  the  effects  of  the 
lifting  surface,  three-dimensional  cavity,  flow  retardation,  etc.  For  a 
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subcavitating  propeller,  an  airfoil  with  the  same  pressure  distribution  and  the  same 

source  distribution  has  been  used^  instead  of  a  supercavitating  cascade.  When  v„ 

2n 

is  represented  as  the  angle  of  the  two-dimensional  velocity  with  respect  to  the  nose- 

tail  line  and  v^n  is  represented  as  the  angle  of  the  three-dimensional  velocity  a 

with  respect  to  the  geometric  advance  angle  3,  then  by  taking  a  single  term  in  the 

summation  of  Equation  (39) ,  aQ  will  be  interpreted  as  the  corrected  angle  of  attack 

of  the  nose-tail  line  with  respect  to  8,  and  c  will  indicate  a  camber  correction  of 

6  ^ 

the  two-dimensional  foil.  In  the  general  case  the  angle  of  attack  will  be 


a 

o 


2 

c 


(40) 


and  the  camber  correction  due  to 


°2J  i xi 


has  to  be  considered  in  addition  to  the 


correction  due  to  cq.  The  optimum  number  of  terms  in  Equation  (39)  depends  upon  the 
number  of  collocation  terms. 


FORCES  ON  THE  BLADE 

As  in  the  two-dimensional  supercavitating  flow,  lift  and  drag  should  be  evalu¬ 
ated  by  integrating  pressure  on  the  blade  surface.  The  section  drag  in  the  direction 
of  the  nose-tail  line  is 

CD  =  c  J  ^p+0''  v3n  dx  +  2  CDV  =  CDC  +  T  CDV 
o 


since 


* 


3n 


dx 


0 


(42) 


6 

where  is  the  friction  drag,  which  is  the  same  as  for  the  subcavitating  case. 
From  the  Kutta-Joukowsky  theorem,  Equation  (22), 


p  +  a 
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V 

s 


V 
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where 


can  now  be  represented  by  induced  velocities. 


Thus 
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By  writing 
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(43) 
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where 
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Likewise  the  sectional  lift  coefficient  can  be  written 


C  =  c 


(p+O)  dx  =  — 
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Since  V£n  and  y  are  available  from  the  preliminary  design  calculations 


Pk  = 


k-1  , 

x  y  dx 


and 


\  = 


k-1  , 

X  V2n  Y  dx 


can  also  be  precalculated  in  the  preliminary  design.  Then  the  nondimensional  thrust 
and  power  coefficients  are  obtained  in  a  conventional  way: 


C„,  = 


T 


T  "  ir  2  2 

fpVg“R 


(C*L  cos  <j)-Cp  sin  <j>)  dr 


C„  - 


(oQ 
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P  ‘  f  p  R2  V3 

2  1  s  r 


LH 
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r(CL  sin  <f>+CD  cos  <f>)  dr 


II 


where  T  =  propeller  thrust 
Q  =  propeller  torque 
4>  =  pitch  angle,,  i.e., 


<|>  =  3  +  a 


(45) 


|«V- 


The  efficiency  of  the  propeller  is  given  by 


C 

n  =  r- 

P 
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NOTE  ON  PRELIMINARY  DESIGN  PROGRAM 

As  for  subcavitating  propellers,  the  preliminary  design  calculations  form  a 
separate  program;  however,  it  is  necessary  to  briefly  explain  this  subject  to  help 
understand  the  final  design  procedures. 

The  present  preliminary  design  of  a  supercavitating  propeller  is  more  or  less  a 
combination  of  subcavitating  lifting-line  theory  and  two-dimensional  supercavitating 
cascade  theory. ^  The  cavity  drag-lift  ratio  will  be  iteratively  fed  from  the  latter 
into  the  former  so  as  to  compute  the  circulation  and  pitch  of  the  blades  needed  to 
produce  a  given  thrust.  In  addition,  the  foil  shape  is  determined  so  that  the  drag- 
lift  ratio  is  not  too  large  for  a  reasonably  thick  cavity.  Both  the  cavity  thick¬ 
ness  near  the  leading  edge,  and  the  minimum  cavity  length  must  be  prescribed  because 
supercavitating  propellers  are  supposed  to  have  a  clean  cavity  covering  the  suction 
side  of  the  blade.  However,  this  cavity  length  is  a  function  of  not  only  foil  shape 
but  also  the  cavitation  number,  which  is  determined  by  the  design  conditions.  To 
design  such  a  foil,  it  is  convenient  to  consider  three  elementary  foils:  a  basic, 
cambered  low-drag  foil,  such  as  the  two-term  camber  foil;  a  flat  plate  to  supply 

angle  of  attack;  and  the  leading-edge  singularity  to  supply  the  leading-edge  cavity 
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thickness.  Called  a  point  drag  because  it  produces  no  lift,  the  leading-edge 
singularity  is  especially  useful  for  creating  a  long  cavity.  This  is  so  because  the 
infinite  cavity  cavitation  number  of  a  supercavitating  cascade  is  linearly  propor¬ 
tional  to  leading-edge  thickness; ^  see  Figure  2.  Figure  2  shows  the  relationship 
between  cavitation  number  and  leading-edge  thickness  for  a  supercavitating  propel¬ 
ler.^  Model  3770  in  Figure  2  will  be  used  for  the  numerical  test  and  is  discussed 
in  more  detail  1.  ter. 

In  general,  the  specification  for  the  leading-edge  cavity  thickness  is  based 
upon  the  strength  requirement.  When  the  cavitation  number  is  relatively  large, 
however,  the  leading-edge  cavity  thickness  may  have  to  be  greater  than  the  strength 
condition  requires.  This  greater  thickness  is  needed  to  allow  cavity  lengths  that 
are  more  than  50  percent  longer  than  the  chord.  Each  section  has  different  cascade 
parameters,  i.e.,  the  solidity  and  the  stagger  angle;^  see  Figure  3.  The  present 
preliminary  design  method  allows  the  option  of  using  either  a  specified  basic  load 
distribution  or  a  specified  basic  camber  shape  for  the  blade  section.  In  an  actual 
design,  the  given  basic  camber  shape  of  the  blade  section  seems  to  be  more  conve¬ 
nient.  Given  the  basic  camber  shape  (the  two-term  camber  in  an  infinite  medium 
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normalized  by  the  lift  coefficient,  for  example)  one  finds  the  shock-free  angle  at 
each  station  of  a  supercavitating  propeller  blade  by  considering  a  supercavitating 
cascade  of  infinite  cavity  length  at  each  blade  section.  The  angle  of  attack  and 
the  point  drag  are  combined  with  the  basic  shock-free  camber  to  supply  the  given 
leading-edge  thickness,  the  minimum  cavity  length,  and  the  required  load  distribu¬ 
tion.  In  this  process  there  are  two  options:  either  the  camber  or  the  angle  of 

11 

attack  per  unit  lift  coefficient  is  preset,  where  the  former  is  Case  1  and  the 
latter  is  Case  2. 

The  hydrodynamic  advance  angle  0^  can  be  determined  by  either  of  two  ways:  any 

form  of  tan  0.  may  be  preset,  such  as  tan  0.  =  c/r,  or  the  optimum  pitch  condition 
l  17  1 

for  the  supercavitating  propeller  may  be  preset 


r-  tan  0 .  = 
A  l 


c(l-t)- 


fa  §)  X 
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where  £  is  the  drag-lift  ratio,  and  t  is  the  thrust  deduction. 

Since  the  local  cavitation  number  is  fixed  according  to  the  design  conditions, 
and  the  lift  coefficient  varies  for  each  iteration,  the  cavity  length  may  also  vary. 
Therefore,  the  sectional  supercavitating  cascade  problem  must  be  solved  for  each 
section  and  for  each  iteration.  This  requires  a  great  deal  of  computer  time.  We 
resolved  this  difficulty  by  conveniently  treating  the  cavity  problem  as  a  foil 
having  infinite  cavity  length  with  a  correction  for  the  finite  cavity  effect."^  The 
former  needs  to  be  calculated  only  once  at  each  blade  section,  regardless  of  the 
number  of  iterations.  Only  the  finite  cavity  effect  is  computed  for  each  iteration, 
assuming  the  same  load  distribution  as  for  infinite  cavity  length.  Thus,  the  finite 
cavity  effect  appears  as  a  reduction  of  angle  of  attack  and  camber. 

The  output  of  the  preliminary  design  program  consists  of  the  thrust  and  power 
coefficients,  circulation  and  lift  distribution,  tan  0^,  pitch  distribution, 
efficiency,  etc.  Normal  velocity  distribution  on  the  cavity  and  foil,  foil  cavity 
shape,  load  distribution,  and  other  necessary  data  are  stored  on  a  tape  to  be  fed 
into  the  lifting-surface  design. 
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NUMERICAL  SCHEMES  AND  COMPUTER  PROGRAM 

There  are  several  computer  programs  for  lifting-surface  design  of  subcavitating 
propellers.  Supercavitating  propellers  are  similar  to  wide-  and  thick-bladed  sub¬ 
cavitating  propellers  when  the  cavity  thickness  is  known;  therefore,  it  would  seem 
reasonable  to  use  an  existing  computer  program  for  subcavitating  propellers  rather 
than  to  start  the  entire  complicated  program  from  scratch.  Among  the  available 
programs,  the  recent  program  by  Kerwin  was  chosen  for  two  reasons:  (1)  it  included 
the  effect  of  sources  and  vortices,  and  (2)  it  could  include  the  effects  of  rake  and 
skew  with  variable  A. 

The  main  differences  between  the  programs  for  supercavitating  and  subcavitating 
propellers  are  as  follows.  The  strengths  of  source  distributions  are  not  known  in 
the  supercavitating  case  because  the  cavity  shape  is  not  known  while  the  blade  thick¬ 
ness  is  assumed  to  be  known  for  the  subcavitating  case.  Therefore,  the  program  for 
solving  the  cavity-source  distribution  is  written  according  to  the  present  theory 
explained  in  the  previous  section  and  becomes  the  main  frame  of  the  present  program. 
Routines  in  the  Kerwin  program  are  used  as  much  as  possible.  The  source  is  distrib¬ 
uted  not  only  on  the  blade  but  also  in  the  cavity  wake  according  to  Equations  (31) 
and  (34).  Therefore,  the  number  of  meshes  for  supercavitating  propellers  is  almost 
as  much  as  two  and  a  half  times  the  number  for  the  subcavitating  case.  The  chord- 
wise  load  distribution  is  taken  from  the  supercavitating  cascade  theory  of  pre¬ 
liminary  design  rather  than  from  airfoil  theory.  The  forces  on  the  blade  are  ob¬ 
tained  from  Equations  (41)  and  (45). 

Figure  4  gives  the  outline  of  the  flow  charts  of  the  lifting-surface  program. 

Attention  should  also  be  drawn  to  the  following  points,  which  differ  from  the 
existing  program  for  subcavitating  propellers. 

Because  of  the  singular  behavior  of  the  kernels  of  integrals  appearing  in  the 
induced  velocity  expressions  in  Equations  (17)  and  (18),  all  the  integrations  related 
to  bound  and  trailing  vorticit>,  and  the  source  distribution  are  integrated  analyt¬ 
ically  in  a  small  interval  of  the  radial  direction  on  the  blade  and  cavity  where  the 
collocation  points  are  located;  using  Equation  (35). 

Because  of  the  singularity  along  the  leading  edge,  both  for  source  and  load 
distribution,  the  interval  that  includes  the  leading  edge  is  treated  separately. 

That  is,  the  leading-edge  load  is  the  integrated  load  obtained  when  the  remaining 
load  distribution  is  subtracted  from  the  total  load.  The  leading-edge  source 
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strength  is  represented  by  the  normal  velocities  of  two-dimensional  cascade  results 
solved  in  the  preliminary  design  program;  see  Appendix  B. 

For  convenience,  short  cavity  effects  of  the  supercavitating  cascade  are  com¬ 
puted  in  the  lifting-surface  theory,  when  the  two-dimensional  cascade  is  used  as  an 
approximate  solution  to  be  corrected. 

After  computing  all  the  B  coefficients  of  Equation  (37),  the  simultaneous 
Equations  (38)  are  solved  by  a  subroutine  which  uses  a  Gaussian  elimination  method. 
The  accuracy  of  the  solution  obtained  by  the  least-squares  method  is  approximately 
checked  by  a  comparison  of  the  left-hand  and  right-hand  sides  of  Equation  (28).  Then 
the  total  source  strength  on  the  blade  is  obtained  from  Equation  (31).  The  velocity 
component  due  to  sources  is  obtained  from  all  the  velocity  component  functions 
created  for  A  in  Equation  (36). 

pq 

Because  many  existing  supercavitating-propeller  models  have  chord  lengths 
smaller  near  the  blade  tip  than  near  the  hub,  angular  intervals  near  the  tip  are  too 
small.  Thus,  the  present  program  allows  intervals  of  all  the  even  degrees,  such  as 
2,  4,  etc.  In  addition,  a  1-deg  interval  is  tested  for  a  six-bladed  propeller.  This 
feature  is  supplied  by  using  the  addition  rule  of  trigonometry,  making  use  of  data 
stored  for  cosine  and  sine  functions. 

When  the  rake  and  skew  are  present  in  the  reference  surface,  the  correct  area 
element  Hdpd<f>  in  all  the  area  integrals  has  to  be  taken  into  account  instead  of 


H.,dpd<t>. 

The  collocation  points  can  be  taken  to  be  10  points  on  the  blade  and  5  points 
on  the  cavity  wake.  However,  if  four  points  are  selected  on  the  blade  this  will  be 
exactly  the  same  as  the  collocation  points  for  the  subcavitating  case,  except  for 
the  extra  points  on  the  cavity  wake. 

In  addition  to  other  preliminary  data,  the  and  in  Equation  (44)  are  pre¬ 
pared  in  the  preliminary  design  and  are  conveyed  through  an  input  tape.  As  in  the 
subcavitating  case,  friction  is  taken  into  account,  but  only  on  the  pressure  side. 

The  output  routine  is  created  in  a  format  similar  to  that  of  Kerwin  by  con¬ 
sidering  data  taken  from  the  supercavitating  cascade  of  the  preliminary  design  for 
the  blade-section  and  the  cavity  shapes. 

The  input  to  the  computer  program  includes  the  same  information  as  is  used  for 
the  subcavitating  propeller  design,  such  as  rpm;  ship  speed;  propeller  diameter;  hub 
diameter;  helical  distance  from  an  arbitrarily  fixed  reference  plane,  distance  from 
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the  x-y  plane  to  the  leading  and  trailing  edges  of  the  blade;  inflow  velocity  ratio 
V^VS’  number  of  collocation  points  on  the  blade,  cavity  wake,  etc.  Also,  as 
mentioned  previously,  the  data  tape  created  by  the  preliminary  design  program  for 
all  the  necessary  sectional  data  from  the  supercavitating  cascade  theory  is  fed  into 
the  lifting-surface  design. 

The  output  consists  of  the  components  of  induced  velocities  resulting  from  the 
vortex  and  source  distributions,  correction  factors  for  source  strength  and  camber, 
angle  of  attack  relative  to  3,  pitch  distribution,  thrust  and  torque  coefficients, 
and  efficiency. 


NUMERICAL  EXPERIMENTS  FOR  SOLUTION 

The  cavity  model,  the  solution  for  cavity  source  strength,  the  number  of  inter¬ 
vals  for  vortex  and  source  distributions,  the  number  of  collocation  points,  and  the 
like  cannot  be  determined  theoretically;  instead,  they  must  be  determined  according 
to  the  behavior  of  numerical  output.  The  output,  especially,  should  show  conver¬ 
gence,  which  could  be  built  into  the  program,  if  the  program  is  simple  and  does  not 
require  too  much  time.  However,  a  large  program,  such  as  the  present  one,  which 
requires  considerable  computer  time,  cannot  be  run  for  all  values  of  parameters  so 
as  to  check  convergence  of  the  solution  for  each  design.  Instead,  it  may  be  enough 
to  check  several  aspects  of  parametric  changes  for  a  typical  case  and  to  assume  that 
the  other  cases  will  reasonably  follow  the  typical  case. 

In  the  present  work,  designing  the  Model  4717  supercavitating  propeller,  various 
convergence  checks  have  been  performed,  and  the  results  will  be  shown  in  the  follow¬ 
up.  The  design  conditions-  for  Model  4717C  are  shown  in  Table  1. 

The.  cavity  truncation  locations  were  varied  from  1.5  to  2.8  chord  lengths  to 
check  the  cavity  model.  Figures  5  and  6  show  that  the  truncation  at  2.2  chord 
lengths  and  at  2.5  chord  lengths  produces  almost  the  same  pitch  and  camber  distri¬ 
butions  (differing  less  than  1%);  the  pitch  distribution  is  about  3%  larger  than  at 
1.8  chord  lengths. 

The  influence  of  the  choice  of  angular  interval  on  pitch  and  camber  distri¬ 
bution  is  also  shown  in  Figures  5  and  6.  By  changing  2-deg  intervals  to  1-deg 

intervals,  the  camber  correction  factor  c  increased  about  6%  and  the  pitch  diameter 
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ratio  P/D  increased  about  1%.  Since  the  magnitude  of  camber  is  so  small,  a  change 
in  the  correction  factor  of  10%  ±s  within  the  manufacturing  error.  If  the  blade  tip 
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has  more  than  five  intervals,  2-deg  intervals  are  sufficient;  finer  intervals,  which 
increase  costs  significantly,  appear  to  be  unnecessary. 

To  check  whether  there  were  enough  collocation  points,  we  chose  10  points  and 
created  the  foil  shape  by  plotting  streamlines.  The  results  have  been  compared  with 
the  case  of  four  collocation  points  given  in  Figures  7  and  8.  The  simple  approxi¬ 
mation  obtained  using  Equation  (39)  is  shown  to  express  amazingly  close  agreement 
with  the  plotted  streamlines. 

In  Equation  (28)  is  the  local  mean  speed,  which  is  not  known  without  exami¬ 
nation,  so  the  value  from  the  lifting  line  theory  has  been  substituted  for  it.  How¬ 
ever,  a  more  accurate  formulation  would  be  Equation  (53)  (Appendix  A)  instead  of 
Equation  (28). 


The  computer  results  for  these  two  cases  were  almost  the  same.  The  present  program 

used  Equation  (47),  although  slightly  more  computer  time  was  needed.  When  G/V  is 
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quite  large,  it  may  improve  the  solution. 

To  check  whether  the  solution  satisfies  the  boundary  conditions  well,  the  left- 
hand  and  right-hand  sides  of  Equation  (47)  were  plotted  in  Figures  9-12;  and  the 
radial  components  of  velocities  on  the  blades  were  plotted,  Figures  13  and  14.  These 
calculations  were  made  with  and  without  the  hub  boundary  conditions  being  satisfied 
when  the  degrees  of  p  and  x  in  the  double  polynomials  in  Equations  (31)  through  (33) 
were  taken  as  3  in  one  case  and  as  5  in  the  other  case.  During  this  process,  we 
noticed  an  interesting  phenomenon:  an  instability  occurred  in  the  numerical  value 
of  radial  velocity  for  the  solution  which  did  not  satisfy  the  hub  boundary  condition. 
That  is,  if  the  hub  boundary  condition  was  not  specified,  a  slight  change  of  parame¬ 
ters,  such  as  cavity  length,  number  of  intervals,  or  the  degree  of  polynomials, 
produced  large  changes  in  radial  velocities.  Yet,  the  numerical  values  of  the  thrust 
and  torque  coefficients  or  the  pitch  distribution  did  not  change  too  much.  This  may 
be  because  the  linear  boundary  conditions  on  a  cavity  or  foil  do  not  include  any 
constraint  on  the  radial  velocity.  In  the  present  problem,  the  only  constraint  on 
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the  radial  velocity  is  on  the  hub  condition.  Thus,  the  hub  boundary  condition  is 
needed  not  only  to  find  the  hub  effect  on  the  pitch  distribution  but  also  to  make 
the  solution  stable. 

The  specified  boundary  conditions  are  well  satisfied  in  general,  although  when 
the  boundary  condition  on  the  hub  is  included,  the  cavity  conditions  are  slightly 
less  accurate,  as  shown  in  Figures  11  and  12.  The  differences  in  the  radial  veloc¬ 
ities  occurring  for  the  case  with  and  the  case  without  the  hub  boundary  condition 
(see  Figures  13  and  14)  indicate  an  instability.  The  radial  velocity  satisfying  the 
hub  boundary  condition  in  Figure  13  is  the  stable  solution.  In  Figure  15  the  pitch- 
diameter  ratio  is  shown  for  Model  4717C  with  hub  images.  In  Figure  16  the  camber 
correction  factors  are  shown  for  Model  4717C  with  and  without  hub  images.  When  the 
computed  results  obtained  without  consideration  of  the  hub  boundary  condition 
happened  to  have  radial  velocities  with  small  values,  the  results  were  very  close  to 
the  solutions  obtained  when  the  hub  boundary  condition  was  considered.  The  numerical 
results  reported  in  the  following  discussion  were  obtained  without  satisfying  the 
hub  boundary  condition.  In  the  cases  given,  the  instability  phenomenon  was  not 
noticed. 


NUMERICAL  EXAMPLES  FOR  PROPELLER  DESIGN  AND  DISCUSSIONS 

Many  supercavitacing  propellers  have  been  designed  and  tested  in  the  past;  of 
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that  number,  two  propellers,  DTNSRDC  Models  3770  and  3870  were  chosen  to  determine 
if  this  design  program  is  reasonable.  The  former  propeller  has  three  blades  and  a 
low  advance  coefficient,  and  the  latter  has  four  blades  and  a  high  advance  coeffi¬ 
cient.  E'.periments  showed  that  both  propellers  had  smooth  cavities.  The  experi- 
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mental  results  and  the  previous  design  calculations  are  available. 

The  design  and  performance  characteristics  of  the  two  propellers  are  shown  in 
Table  2. 

It  is  extremely  difficult  to  compare  the  present  numerical  results  to  the 
experimental  results  for  propellers  that  were  designed  using  an  entirely  different 
method.  The  present  program  is  intended  for  design,  not  prediction.  The  present 
program  does  not  produce  data  on  leading-edge  cavity  thicknesses,  input  that  is 
essential  to  the  design  of  propellers  similar  to  Models  3770  and  3870.  To  check 
the  reasonableness  of  the  present  program  we  guessed  at  the  leading-edge  cavity 
thicknesses  for  Models  3770  and  3870;  this  is  presented  in  Figure  17.  Because  the 
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actual  cavity  thickness  was  never  measured,  this  is  not  a  scientific  estimate.  The 
leading-edge  cavity  thickness  selected  for  Model  3770  results  in  an  almost  infinite 
cavity  length  at  every  section  of  the  blade,  except  near  the  tip  and  the  hub.  How¬ 
ever,  because  the  design  cavitation  number  of  Model  3870  is  not  very  small,  the 
selected  leading-edge  cavity  thickness  is  not  large  enough  to  induce  a  smooth  sheet 
cavity  all  over  the  blade.  Therefore,  the  leading-edge  cavity  thickness  of  Model 
3870  was  corrected  to  give  a  cavity  length  at  least  50  percent  longer  than  the  chord. 
To  do  this,  an  extra  leading-edge  point  drag  was  added  to  the  cascade  theory  used  in 
the  present  design  method,  as  explained  previously.  The  design  thrusts  used  in  the 
present  lifting-line  computations  are  the  experimental  values  listed  in  Table  2, 
where  Case  1  is  for  cases  without  angle  of  attack  and  Case  2,  with  preset  angle  of 
attack.  It  is  seen  that  the  efficiencies  of  the  propellers,  as  predicted  by  the 
present  method,  are  very  close  to  the  measured  efficiencies,  even  in  the  preliminary 
design  stage  of  calculations. 

Pitch  distributions  obtained  from  the  preliminary  design  and  lifting-surface 
design  computations,  according  to  the  two  design  approaches,  Case  1  (without  angle 
of  attack)  and  Case  2  (with  angle  of  attack)  are  shown  in  Figures  18  through  21 
together  with  the  pitches  of  the  two  propeller  models.  The  pitch  values  obtained 
from  the  preliminary  design  calculations  are  higher  than  those  obtained  from  lifting- 
surface  calculations  because,  in  preliminary  design,  the  effect  of  flow  retardation 
is  not  considered.  The  pitch  distribution  is  also  related  to  the  leading-edge 
cavity  thickness.  In  general,  when  the  leading-edge  thickness  increases,  the  pitch 
also  increases;  however,  the  efficiency  decreases  slightly.  The  pitch  distributions 
for  the  predictions  and  for  the  models  are  noticeably  different.  This  is  because 
the  optimum  lift  distribution  and  pitch  angle,  which  are  influenced  by  the  blade- 
cavity  interference  in  the  present  method  are  quite  different  from  those  of  the 
models.  If  these  factors  are  taken  into  account,  all  results  appear  reasonable 
compared  with  those  of  the  models. 

The  lifting-surface  corrections  to  the  source  distribution  for  the  two  pro¬ 
pellers  ar<  shown  in  Figure  22.  Although  the  correction  for  Model  3770  is  close  to 
1,  it  is  about  20  percent  greater  near  the  trailing  edge  than  at  the  leading  edge  for 
Model  3870.  If  the  correction  factor  is  unity,  this  means  that  the  cascade  source 
strength  is  the  same  as  the  blade  cavity  source  strength.  This  source  strength  has 
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the  main  influence  on  pitch,  thrust,  and  efficiency.  The  calculated  efficiency  is 

close  to  that  obtained  in  the  model  experiments,  as  shown  in  Table  2. 

The  lifting-surface  camber  correction  factors  for  the  two  propellers- are  shown 

in  Figures  23  and  24.  The  camber  correction  factor  for  Model  3870  is  much  larger 

2 

than  for  Model  3770,  as  in  the  calculation  of  Venning  and  Haberman.  The  pitch 
distribution  and  camber  correction  curves  for  Model  3870  are  quite  different  from 
those  for  Model  3770.  The  former  has  shorter  cavities  with  four  blades  and  large 
expanded  area  ratio  (EAR),  and  the  latter  has  longer  cavities  with  three  blades  and 
smaller  EAR. 


EXPERIMENTAL  EVALUATION  OF  THEORY 

Finally,  the  Center  conducted  an  experimental  program  to  evaluate  the  present 

method  for  designing  supercavitating  propellers.  Two  supercavitating  propellers 
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were  designed,  using  the  present  method,  for  a  200- ton  (181  metric  ton)  hydrofoil 
craft.  The  propeller  design  criteria  are  given  in  Table  3.  The  propeller  design 
characteristics  are  given  in  Table  2. 

Two  model  propellers  were  manufactured  from  these  designs.  The  geometry  of 
these  propellers  .is  given  in  Table  4  and  drawings  of  the  propellers  are  shown  in 
Figures  25  and  26. 

The  experimental  program  was  divided  into  two  phases.  The  results  of  the  first 
phase,  the  measurements  of  blade-cavity  shapes,  are  reported  in  the  following 
section,  while  the  results  of  model  propeller  performance  are  described  in  the  last 
section. 


EXPERIMENTAL  INVESTIGATION  OF  BLADE-CAVITY 
THICKNESS  DISTRIBUTION 

A  series  of  experiments  was  performed  to  determine  how  well  linear  theory 

predicted  the  upper  cavity  surface  location  for  Propellers  47(-7C  and  4738A.  For 

these  experiments,  brass  pins  of  varying  lengths  were  attached  to  the  backs  of  the 

propeller  blades.  During  propeller  operation  in  the  36-in.  variable-pressure  water 

tunnel,  one  could  see  when  the  pins  came  into  contact  with  the  upper  cavity  surface. 

This  experimental  procedure  has  already  been  used  to  verify  the  upper  cavity 
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surface  location  for  Propeller  4699,  ’  and  the  parent  design  of  Propellers  4738A 


and  4717C.  In  addition,  cavity  heights  have  been  measured  in  similar  ways  for  a 
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supercavitating  flat  plate  and  a  Tulin-two-term  section  by  Christopher  and  Johnson. 

DESCRIPTION  AND  LOCATION  OF  PINS 

Number-four  brass  machine  screws  with  heads  cut  off  were  used  as  pins.  These 
pins  were  attached  to  the  blade  backs  by  drilling  and  tapping  holes  perpendicular  to 
the  surface  at  12  locations.  The  upper  cavity  surface  location  was  defined  as  a 
point  on  a  line,  perpendicular  to  the  nose  tail  line,  that  runs  through  the  center 
of  the  tapped  hole  at  the  blade  surface;  see  Figure  27.  The  holes  were  drilled  and 
tapped  perpendicular  to  the  back  of  the  blade  to  cant  the  pins  slightly  away  from  the 
reference  line  (the  line  perpendicular  to  the  nose  tail  line).  A  slight  error  was 
introduced,  but  the  machining  process  was  greatly  simplified.  The  locations  of  these 
pins  were  as  follows:  10,  30,  60,  and  90  percent  of  chord  at  nondimensional  radii 
(r/R)  of  0.361,  0.544,  and  0.726;  see  Figure  28.  Since  the  blade  at  10  percent  of 
chord  was  too  thin  to  tap,  the  pins  at  these  locations  were  soldered  in  place.  The 
pins  at  all  other  locations  were  screwed  in  and  secured  by  tiny  electrical  lock  nuts 
(see  photos  in  Reference  19) . 

Three  sets  of  pins  were  used  for  testing  Propeller  4717C.  When  installed,  the 
first  set  of  pins  protruded  above  the  back  of  the  blades  to  a  height  that  corre¬ 
sponded  to  three  times  the  distance  from  the  back  of  the  blades  to  the  theoretically 
predicted  upper  cavity  surface.  The  second  set  of  pins  protruded  by  a  factor  of  1.67 
and  the  third  set  by  a  factor  of  1.0,  the  latter  being  the  theoretically  predicted 
cavity  height.  However,  the  pins  at  10  percent  of  chord  varied  from  this  order; 
their  heights  corresponded  to  cavity-height  factors  of  3,  2.33,  and  1.0.  The  pins  at 
the  10  percent  of  chord  locations  had  to  be  filed  by  hand'  to  the  correct  height,  and 
the  factor  of  2.33  rather  than  1.67  was  used  to  ensure  that  the  pin-height  would 
exceed  the  experimental  cavity  thickness.  The  experiments  have  shown,  however,  that 
the  theory  overpredicted  the  cavity  height  near  the  leading  edge,  and  the  pins  could 
have  been  filed  to  a  height  corresponding  to  a  factor  of  1.67. 

Four  sets  of  pins  were  used  for  Propeller  4738A.  The  first  set  of  pins  pro¬ 
truded  above  the  back  of  the  blade  to  a  height  corresponding  to  a  factor  of  1.8  times 
the  distance  from  the  back  of  the  blade  to  the  theoretically  predicted  cavity  sur¬ 
face.  The  other  three  sets  corresponded  to  factor  of  1.4,  1.0  (theoretical)  and 


0.6.  The  pins  at  the  10  percent  of  chord  locations  were  filed  to  a  height  corre¬ 
sponding  to  the  multiplication  factor  used  for  the  other  pins  in  the  set  being 
tested. 

A  thin  coat  of  international  yellow  paint  was  applied  to  the  tips  of  the  pins 
prior  to  testing  to  aid  in  visual  observation.  The  experimenters  could  thus  locate 
the  pins  easily  when  the  propeller  was  revolving. 

EXPERIMENTAL  PROCEDURE 

The  propeller  rpm  was  increased  until  the  cavity  enclosed  all  the  pins.  This 
procedure  was  begun  with  the  '  t  of  brass  pins  that  protruded  highest.  While  the 
pressure  and  velocity  in  the  36-in.  variable  pressure  water  tunnel  was  held  constant, 
corresponding  to  the  design  cavitation  number,  a  =  0.34,  the  propeller  rpm  was 
gradually  reduced  in  decrements  of  10.  Each  time  the  propeller  rpm  was  reduced,  a 
hand-held  strobe  unit  was  used  to  observe  visually  whether  the  pins  were  in  or  out 
of  the  cavity.  These  observations  could  be  made  rapidly,  although  the  pins  that 
barely  touched  the  upper  cavity  surface  required  more  attention  than  did  the  others. 
When  all  pins  protruded  through  the  cavity  surface,  the  procedure  was  repeated  with 
another  set  of  pins.  A  depth  micrometer  was  used  before  and  after  each  test  to 
measure  the  height  of  the  pins  above  the  back  of  each  blade.  This  was  done  to  ensure 
that  the  pins  were  at  the  correct  height  and  had  not  moved  during  testing. 

As  the  pins  began  to  break  through  the  cavity  surface,  a  furrow  or  small 
groove  formed  in  the  surface,  accompanied  by  some  spray  or  cavitation  behind  the 
pin.  A  directional  strobe  unit  with  variable  light  intensity  made  these  furrows 
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much  more  visible.  This  phenomenon  has  been  recorded  in  several  color  photographs. 

Large,  international-yellow  numbers  painted  on  the  backs  of  the  blades  proved 
invaluable  during  testing.  Also,  each  propeller  hub  was  coded  with  a  series  of  dots, 
the  number  of  which  corresponded  to  the  number  on  each  blade. 

DISCUSSION  OF  EXPERIMENTAL  RESULTS 

Figures  29-31  compare  linear  theory  predictions  of  cavity  height  with  experi¬ 
mental  results  for  Propeller  4717C.  These  figures  show  three  experimental  upper 
cavity  surfaces  corresponding  to  three  values  of  J,  one  of  which  is  the  design  J 
(1.037).  Note  that  the  following  relationship  gives  the  advance  angle,  3,  at  each 
radial  blade  section, 
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(VA/2irrn) 


3  =  tan 

=  tan  1  [J/(rrr/R)] 

From  this  relationship,  one  can  determine,  approximately,  the  corresponding  shifts 
in  upper  cavity  surface  that  result  from  small  changes  in  angle  of  attack.  For 
example,  in  Figure  29,  at  the  design  J  (1.037),  3  is  determined  to  be  42.44  deg. 

For  a  J  value  of  1.0,  3  -  41.40  deg.  Therefore,  a  change  of  J  corresponding  to 
0.037  has  caused,  approximately,  a  one-deg  change  in  angle  of  attack,  which  shifts 
the  cavity  surface  upwards  as  shown. 

Theoretically,  the  section  lift  and  cavity  thickness  for  Propeller  4717C  are 
generated  entirely  by  camber  and  point  drag  (note  the  blunt  nose  in  Figures  29 
through  31).  That  is,  no  incidence  was  used  in  the  design  to  generate  lift  or  cavity 
thickness.  In  Figures  29  though  31,  the  theoretical  prediction  of  cavity  height 
agrees  fairly  well  with  the  experimental  data.  Near  the  leading  edge,  however,  the 
theory  appears  to  overpredict  cavity  thickness.  Also,  visual  observations  indicated 
that  the  backs  of  the  blades  at  all  radial  sections  on  Propeller  4717C  were  wetted 
to  about  2-  or  3-percent  of  chord  from  the  leading  edge.  At  this  point,  separation 
was  caused  by  a  locally  flat  area  that  was  inadvertently  machined  onto  the  back  of 
the  blade.  Although  this  local  flat  was  almost  microscopic,  it  effectively  caused 
separation.  Apparently,  very  near  the  leading  edge  some  portion  of  the  blade  metal 
was  interfering  with  the  upper  cavity  streamline. 

Figures  32  to  34  compare  linear  theory  predictions  of  cavity  height  with  experi¬ 
mental  results  for  Propeller  4738A.  Note  in  these  figures  the  large  amount  of  point 
drag  or  blunt  nose  indicated  by  the  theory.  This  results  because  both  Models  4738A 
and  4717C  were  designed  to  have  approximately  the  same  full-scale  stress  levels. 

This  dictated  that  the  maximum,  theoretical,  cavity  thicknesses  for  Models  4738A  and 
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4717C  would  be  almost  the  same.  To  obtain  the  same  maximum  thickness  in  a  shorter 
distance,  we  used  a  large  amount  of  point  drag  together  with  incidence  and  camber  to 
generate  the  theoretical  cavity. 

Figures  32  to  34  show  three  experimental  upper  cavity  surfaces  corresponding  to 
three  values  of  J.  Note  that  at  r/R  =  0.361,  Figure  32,  the  blade  was  fully  wetted 
at  the  design  value  of  J  (1.037);  therefore,  cavity  heights  for  three  other  values 
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of  J  have  been  shown.  At  a  J  value  of  0.98,  the  experimental  cavity  surface  coin¬ 
cides  with  the  predicted  cavity  surface  height  at  aft  locations  on  the  blade.  This 
J  value  represents  an  approximate  increase  in  angle  of  attack  of  1.6  deg  over  the 
design  value  of  J.  For  the  cavity  surface,  corresponding  to  a  J  value  of  0.96,  the 
increased  incidence  is  equal  to  about  2.2  deg.  However,  as  mentioned  before,  the 
cavity  did  not  spring  from  the  leading  edge.  It  moved  down  the  span  as  rpm  was  in¬ 
creased.  Since  the  water-tunnel  velocity  was  held  close  to  35  fps  (10.688  m/sec) 
and  since  the  model  propeller  diameter  was  16  in.  (40.64  cm),  the  difference  in  model 
rpm  corresponding  to  the  J  values  of  0.98  and  1.037  was  88.  This  corresponds  to  an 
increase  of  about  58  in  full-scale  rpm.  It  is  also  interesting  to  note  that,  accord¬ 
ing  to  performance  evaluation  experiments,  an  increase  of  29  rpm,  over  the  1000  rpm 
of  full-scale  design  would  give  the  design  thrust. 

At  the  two  outer  radial  positions,  r/R  =  0.544  and  0.726,  full  cavitation  did 
occur  at  the  design  J  (1.037),  but  the  theory  overpredicted  the  cavity  surface 
height.  As  with  Propeller  4717C,  the  back  of  the  blade  near  the  leading  edge  of 
Propeller  4738A  was  wetted  to  about  2-  or  3-percent  of  chord. 

To  understand  more  fully  the  discrepancy  between  theory  and  experiment,  the 
reader  should  recall  that  a  point  drag  is  a  linear  theoretical  model  of  the  leading 
edge  cavity  thickness  represented  by  a  point  singularity.  Experimental  results  indi¬ 
cate  that  the  actual  separation  point  at  the  leading  edge  must  be  carefully  chosen, 
for  example,  as  a  slope  discontinuity  of  the  blade  surface,  to  achieve  the  designed 
leading  edge  cavity  thickness;  if  the  predicted  leading  edge  cavity,  not  just  70 
percent  of  it,  had  been  filled  by  a  material  up  to  2  percent  of  chord  from  the 
leading  edge,  the  experimental  results  would  have  almost  coincided  with  the  theory, 
except  very  near  the  hub,  where  the  hub  effect  is  important. 

CAVITATION  PERFORMANCE  CHARACTERISTICS  OF  SUPERCAVITATING 
PROPELLERS  4717B,  4717C  AND  4738A 

BACKGROUND 

Propeller  4717C  was  originally  manufactured  as  Propeller  4717B.  Propeller  4717B 
was  identical  to  Propeller  4717C  except  for  the  backs  of  the  blades,  which  had  a 
shape  to  conform  to  the  predicted  cavity  shape  at  design  operating  conditions.  The 
primary  purpose  of  Propeller  4717B  was  to  determine,  by  observation,  how  well  the 


blade  section  shape  represented  the  blade  cavity  shape.  Following  characterization 
and  observation,  Propeller  4717B  was  finish  cut  to  the  final  design  version,  Pro¬ 
peller  4717C. 

Propeller  4738A  is  a  six  bladed  propeller  with  the  same  expanded  area  as  the 
previous  four  bladed  propeller,  and  is  designed  for  the  same  conditions  as  the  other 
propellers. 

EXPERIMENTAL  PROCEDURE 

Cavitation  performance  characteristics  and  cavitation  observations  were  ob¬ 
tained  in  the  36-in.  variable  pressure  water  tunnel.  Tunnel  water  velocities  were 
measured  by  the  tunnel  venturi  system.  The  scope  of  the  experiments  is  given  in 
Table  5. 

Tunnel  pressure  and  water  velocity  were  set  to  establish  each  cavitation  number 
and  then  propeller  revolution  rate  was  varied  to  cover  a  range  of  advance  coeffi¬ 
cients.  Propeller  thrus  and  torque  were  measured  at  each  condition  and  sketches 
were  made  of  the  cavitat,  n  present.  The  Reynolds  number,  R  ,  during  the  experiments 

C  /■  ft 

ranged  from  7.5  x  10  to  5.6  x  10  . 

PRESENTATION  OF  DATA  AND  DISCUSSION 

The  thrust  and  torque  data  were  reduced  to  nondimensional  coefficients  of  KT  and 
Kq.  Propeller  efficiencies  were  calculated  from  faired  values  of  K^,  and  K^.  The 
cavitation  performance  characteristics  of  the  three  propellers  are  presented  in 
Tables  C  through  8. 

Curves  representing  the  faired  data,  from  Tables  6  through  8,  are  shown  as  an 
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example  in  Figure  35.  Curves  of  maximum-speed  thrust  loading  (K  /J  )  have  been 
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added  to  the  performance  curves  for  Propeller  4738A.  The  intersection  of  the  K^J 
curve  and  the  K^,  curves  at  the  design  sigma  (a)  determines  the  predicted  operational 
point  for  each  propeller.  A  comparison  between  the  design  operational  points  and 
the  points  predicted  by  the  experimental  data  is  given  in  Table  9. 

Sketches  of  the  back  cavitation  present  on  the  propellers  at  two  cavitation 
numbers  are  given  in  Figures  36  through  38.  These  sketches  cover  a  range  of  advance 
coefficients  from  partially  cavitating  to  fully  cavitating  conditions.  If  advance 
coefficients  had  lower  values  than  those  shown,  the  propellers,  at  the  same  cavi¬ 
tation  number,  would  also  be  fully  cavitating.  With  only  one  exception,  propellers 
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contained  no  face  cavitation  over  the  range  of  cavitation  numbers  and  advance  co¬ 
efficients  covered.  The  sole  exception  was  at  an  advance  coefficient  of  1.2,  where 
some  leading  edge  face  cavitation  was  observed  at  cavitation  numbers  of  0.75  and 
lower . 


CONCLUSIONS 

At  design  speed  coefficient  and  design  a,  Propeller  4717B  contains  practically 
no  back  cavitation.  If  advance  coefficient  is  reduced  slightly,  at  design  a,  the 
backs  of  the  blades  are  covered  by  sheet  cavitation  from  the  blade  tip  to  50  percent 
radius.  This  indicates  that  the  predicted  cavity  shape  over  this  part  of  the  pro¬ 
peller  blades  is  quite  accurate. 

Neither  of  the  designed  propellers,  4717C  and  4738A,  had  face  cavitation  at  the 
design  operational  points.  Propeller  4717C  essentially  had  full  back  cavitation  and 
Propeller  4738A  had  back  cavitation  from  about  35  percent  radius  to  the  tip  of  the 
blades  at  the  design  operational  point. 

The  propeller  theory  slightly  overpredicted  the  available  thrust  for  both  pro¬ 
pellers.  Propeller  4717C  would  require  6  percent  more  rpm  and  8  percent  more  power 
than  predicted  to  reach  design  speed.  Propeller  4738A  would  require  5  percent  more 
rpm  but  4  percent  less  power  than  predicted  to  reach  design  speed.  If  the  operating 
point  is  defined  as  the  speed  and  rpm  where  the  propellers  absorb  the  available 
maximum  power,  Propeller  4717C  would  operate  at  =  58.7  knots  and  rpm  =  1016,  and 
Propeller  4738A  would  operate  at  =  61  knots  and  rpm  =  1054.  The  propeller 

efficiencies  at  these  conditions  are  66  percent  and  67  percent,  respectively. 
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It  has  been  recognized  that  the  nonlinear  effects  on  lift  and  drag  of  cavi- 

2 

tading  foils  are  approximately  equal  to  -0.5  C^/(l+o)  and  -0.5  CpC^/(l+a).  There¬ 
fore,  propellers  designed  according  to  the  linear  theory  would  produce  less  thrust, 
as  indicated  in  the  experimental  results.  However,  the  leading  edge  cavity  thickness 
was  slightly  smaller  than  the  design  thickness  due  to  the  unmatched  separation  point 
and  70  percent  filling  of  cavity  thickness.  Thus  the  drag  should  have  been  a  little 
less  than  predicted  by  the  linear  theory,  as  was  found  experimentally  for  Propeller 
4738A.  Because  of  the  decreased  efficiency  of  Propeller  4717C,  however,  more  care 
may  be  required  in  calculating  the  cavity  drag  due  to  the  blunt  leading  edge. 

All  in  all,  the  design  theory  predicted  the  cavity  thickness  and  the  propeller 
performance  quite  well  at  the  design  point,  within  the  bounds  of  error  to  be 
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expected  of  linear  theory.  It  should  be  remembered  that  the  present  design  theory 
does  not  include  semiempirical  formulae  or  factors.  In  addition,  this  verifies  that 
the  point  drag  of  linear  theory  is  very  useful  in  solving  one  of  the  important 
problems  relating  to  supercavitating  propellers — how  to  make  propeller  blades  having 
sufficiently  thick  leading  edges  without  paying  too  much  of  an  efficiency  penalty. 
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(DESIGN  SECTION 
ff/CL) 

r/R  =  0.25  J 


ul 


0. 5/5.6^ 


LEADING  EDGE  THICKNESS,  t/CL  c(x.,  =  0.1) 

Figure  2  -  Example  of  Relation  Between  Choked  Flow 
Cavitation  Number  and  Leading-Edge  Thickness  of 
Supercavitating  Cascade  for  Blade  Sections  of 
Supercavitating  Propeller  Model  3770 


RADIUS,  (r/R) 

Figure  3  -  Example  of  Cascade  Parameters  for 
Supercavitating  Propeller  Model  3770 
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Figure  4  -  Flow  Chart  for  Lifting-Surface  Design  of 
Supercavitating  Propellers 
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/✓ 

Figure  5  -  Computed  Camoer  Correction  Factors  cq  of  Supercavitatmg 

Propeller  Model  4717  for  Different  Angular  Intervals  <5  and 
Cavity  Truncation  Points  £. 


Figure  6  -  Computed  Pitch-Diameter  Ratio  of  Supercavitating  Propeller 
Model  4717  for  Different  Angular  Intervals  6  and  Cavity  Truncation 

Points  Z 


0.2  0.4  0.6  0.8  1.0 

x/c 

Figure  7  -  Computed  Shape  of  Blade  Face  of  Model  4717  at  r/R  =0.8 
for  Two  Numbers  of  Collocation  Points 
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WITH  EQUATION  (39) 


Figure  8  -  Computed  Shape  of  Blade  Face  of  Model  4717  at  r/R  =  0.3978 
for  Two  Numbers  of  Collocation  Point 
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MODEL  3770 


r/R  =  0.3 


MODEL  3870 


Figure  26  -  Drawing  of  Propeller  4738A 


Figure  27  -  Placement  of  Brass  Pin  Perpendicular  to  Back  of  Blade 
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TABLE  1  -  PROPELLER  DESIGN  PREDICTIONS 


Propeller  Number 

4717C 

4738A 

V 

s 

58.84  ’ 

58.84  knots  (30.35  m/sec) 

n 

1000 

1000 

*T 

0.092 

0.092 

J 

1.037 

1.037 

Z 

4 

6 

P/D0.7 

1.416 

1.410 

0 

0.34 

0.34 

n 

0.67 

0.67 

TABLE  2  -  DESIGN  AND  PERFORMANCE  CHARACTERISTICS  OF 
SUPERCAVITATING  PROPELLERS 


TABLE  3  -  PROPELLER  DESIGN  CRITERIA 


Power,  maximum  continuous 

16,000  hp  ('11,931  kW) 
for  two  propellers 

1  -  t 

0.925 

1  -  wt 

0.870 

Shaft  Centerline  at 

Design  Speed 

6.82  ft  (2.079  m) 

Maximum  D 

5  ft  (1.524  m) 

No.  of  Blades 

4,  6 

Design  Objective 

maximum  speed 

Minimum  Rotative  Speed 

750  rpm 

Hump  Thrust  Margin 

20  percent 

TABLE  4  -  MODEL  PROPELLER  GEOMETRY* 


Design  Parameters 

Propeller  Number 

4717C 

4738A 

Z 

4 

6 

Diameter,  in.  (cm) 

16.000  (40.64) 

16.000  (40.64) 

Pitch,**  in.  (cm) 

22.656  (57.546) 

22.560  (57.302) 

P/D** 

1.416 

1.410 

EAR 

0.495 

0.492 

Chord  length,**  in.  (cm) 

4.560  (11.582) 

3.017  (7.663) 

interpreted  from  drafting  room  offsets. 

**At  0.7R. 

TABLE  5  -  SCOPE  OF  EXPERIMENTS 


Propeller  Number 

4717B 

4717C 

4738A 

0 

V* 

J 

0 

V* 

J 

a 

V* 

J 

6.0 

15.0 

0.45-1.3 

6.0 

15.0 

0.45-1.3 

3.0 

30.0 

0.5-1. 2 

3.0 

30.0 

0.5-1. 2 

3.0 

30.0 

0.5-1. 2 

1.5 

30.0 

0.5-1. 2 

1.5 

30.0 

0.5-1. 2 

1.5 

30.0 

0.5-1. 2 

0.75 

30.0 

0.5-1. 2 

0.75 

30.0 

0.5-1. 2 

0.75 

30.0 

0.5-1. 2 

0.5 

30.0 

0.6-1. 2 

0.3 

35.0 

0.7-1. 2 

0.34 

35.0 

0.7-1. 2 

0.34 

35.0 

mu 

*In  feet  per  second. 
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TABLE  6  -  CAVITATION  PERFORMANCE  OF  PROPELLER  4717B 


SIGMA  (a)*  =  3.000 


SIGMA*  *  1.500 


J 

Kt 

10KQ 

n 

J 

kt 

10KQ 

n 

.5000 

.2014 

.4583 

.3497 

1 

.2808 

.3213 

.5500 

.2213 

.5185 

.3744 

■ 

.2997 

.3576 

.6000 

.  2329 

.5450 

.4080 

.1421 

.3388 

.4003 

.  6500 

.2334 

.5491 

.4398 

.6500 

.1623 

.3802 

.4417 

.7000 

.  2238 

.5395 

.4622 

.  7000 

,1770 

.4126 

.4778 

.7500 

.2055 

.5219 

.4701 

.  7500 

.1826 

.4302 

.5068 

.6000 

.  1907 

.5000 

.4602 

.6000 

.1783 

.4312 

.5254 

.  0.600 

.  1520 

.4755 

.4324 

.0500 

.1648 

.4174 

.5339 

.0000 

.  1  220 

.  449 1 

.3893 

.9000 

.1440 

.3925 

.5256 

.0500 

.0931 

.4203 

.3350 

.9500 

.1135 

.3608 

.4969 

1 .0000 

.0670 

.  3885 

2743 

1 .0900 

.0910 

.3267 

.4434 

t  .0500 

.0442 

.3533 

.2090 

1  .0500 

.0633 

.2931 

.3608 

1 . tooo 

.0241 

.3146 

.1341 

1.1000 

.0360 

.2602 

.2424 

1.1500 

.0042 

.2736 

.0263 

1.1500 

.0082 

.2249 

.0671 

♦.0000 

-.0199 

.2329 

-.1632 

1.2000 

-.0236 

.1789 

-.2515 

SIGMA* 

=  0.750 

SIGMA* 

=  0.340 

.5000 

.0920 

.2333 

.3138 

.7000 

.0967 

.2516 

.4281 

.5500 

.0968 

.2525 

.3357 

.7500 

.  1024 

.2663 

.4592 

.6000 

.0987 

.2600 

.3623 

.3000 

.0941 

.2504 

.4786 

.6500 

.1020 

.2668 

.3955 

.8500 

.0854  • 

.2337 

.4943 

.7000 

.'1063 

.2778 

.4344 

.9000 

.0802 

.2261 

.5081 

.7500 

.1168 

.2932 

.4753 

.9500 

.0761 

.2241 

.5134 

.6000 

.1251 

.3102 

.5133 

1 .0000 

.0677 

.2173 

.4957 

.0500 

.1302 

.3241 

.5433 

1 .0500 

.0499 

.1950 

.4273 

.9000 

.1289 

.3596 

.5603 

i’.  1000 

.0214 

.1528 

.2451 

.9500 

.1108 

.3220 

.5579 

1.1500 

-.0119 

.0992 

-.2194 

1 .0000 

.  0987 

.2988 

.5259 

1 .2000 

-.0337 

.0617 

t  .0500 

.0698 

.2610 

.4466 

♦.1000 

.0357 

.2142 

.2920 

1.1500 

.0041 

.1699 

.0439 

1.2000 

-.0134 

.1473 

-.1742 

*Cavitation  number. 
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TABLE  7  -  CAVITATION  PERFORMANCE  OF  PROPELLER  4717C 


SIGMA* 

=  3.000 

SIGMA* 

=  1.500 

J 

kt 

10KQ 

0 

J 

n 

10KQ 

n 

.0000 

.2036 

.4625 

.36o/ 

.5000 

m 

.2813 

.3283 

.5500 

.2201 

.5000 

.3853 

.5500 

.3078 

.3688 

.6000 

.2379 

.5366 

.4218 

.3455 

.4073 

.6500 

.2496 

.5666 

.4562 

.6500 

.3843 

.4453 

.7000 

.2522 

.5780 

.4661 

.7000 

.4171 

.4830 

.7500 

.2443 

.5718 

.5099 

.7500 

HM9 

.4397 

.5199 

.6000 

.2274 

.5504 

5261 

.8000 

.I960 

.4498 

.5547 

.8500 

.2044 

.5182 

.5337 

.8500 

.1935 

.4472 

.5853 

.9000 

.1787 

.4808 

.5325 

.9000 

.1839 

.4328 

.60b7 

.9500 

.1536 

.4433 

.5240 

.9500 

.1678 

.4087 

.6207 

1.0000 

.1314 

.4092 

.5109 

1.0000 

.1461 

.3773 

.6161 

1.0500 

.1126 

.3795 

.4960 

1.0500 

.1204 

.3413 

.5894 

1.1000 

.0955 

.3507 

.4769 

1.1000 

.0927 

.3029 

.5358 

1 . 1500 

.0750 

.3143 

.4369 

1.1500 

.0655 

.2635 

.4549 

1  .2000 

.0420 

.2551 

.3143 

1.2000 

.0416 

.2234 

.3557 

SIGMA* 

*  0.750 

SIGMA* 

=  0.  3t0 

.5000 

.0949 

.2331 

.3239 

.2511 

.4595 

.  5500 

.1056 

.2569 

.3600 

.2488 

.4956 

.6000 

.1059 

.2623 

.3857 

ill  :  c^| 

.241  1 

.5086 

.6500 

.1055 

.2644 

.4129 

.2305 

.5236 

.7000 

.1093 

.271 1 

.4492 

.9000 

.2203 

.5526 

.7500 

.1184 

.2849 

.4960 

.9500 

.0840 

.2137 

.5941 

.8000 

.1311 

.3041 

.5490 

1.0000 

.0843 

.2120 

.6327 

.  CSOO 

.  440 

.3242 

.6009 

1.0500 

.0825 

.2132 

.6469 

.9000 

.1528 

.3394 

.6450 

1.1000 

.0746 

.2108 

.6193 

.9500 

.  1  537 

.3441 

.6755 

1.1500 

.0561 

.1924 

.5335 

1.0000 

.1439 

.3338 

.6861 

1.2000 

.0233 

.1381 

.3220 

1.0500 

.1230 

.3072 

.6691 

1.1000 

.0939 

.2671 

.6156 

1.1500 

.0638 

.2218 

.5262 

1.2000 

.0451 

.1869 

.4609 

*Cavitation  number. 
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TABLE  8  -  CAVITATION  PERFORMANCE  OF  PROPELLER  473 8A 


SIGMA* 

=  6.000 

SIGMA* 

=  3.000 

J 

kt 

10KQ 

0 

J 

kt 

10KQ 

n 

.4500 

.2358 

.5204 

.3245 

.6000 

.2350 

.5097 

.4404 

.5000 

.2808 

.6114 

.3654 

.6500 

.2449 

.5327 

.4755 

.5500 

.3029 

.6592 

.4022 

.7000 

.2530 

.5510 

.5115 

.6000 

.3077 

.6752 

.4352 

.7500 

.2542 

.5583 

.5434 

.6500 

.3000 

.6687 

.4641 

.8000 

.2463 

.5520 

.5681 

.7000 

.2837 

.6472 

.4884 

.8500 

.2297 

.5327 

.5833 

.7500 

.2621 

.6163 

.5077 

.9000 

.2065 

.5032 

.5877 

.8000 

.2377 

.5799 

.5218 

.9500 

.1794 

.4670 

.5807 

.8500 

.2123 

.5408 

.5311 

1.0000 

.1513 

.4278 

.5628 

.9000 

.  1873 

.5005 

.5360 

1’.  0500 

.1242 

.3880 

.5350 

.9500 

.1633 

.4596 

.5371 

1.1000 

.0986 

.3480 

.4960 

1.0000 

.1404 

.4181 

.5344 

1.1500 

.0725 

.3052 

.4346 

1  .0500 

.1183 

.3753 

.5270 

1.2000 

.0405 

.2525 

.3066 

1.1000 

.0963 

.3302 

.5104 

1.1500 

.0730 

.2818 

.4739 

1.2000 

.0468 

.2291 

.3904 

1.2500 

.0159 

.1715 

.1851 

1  3000 

-.0219 

.  1048 

-.4188 

SIGMA* 

=  1.500 

SIGMA* 

=  0.750 

.2497 

.3194 

.5000 

.0841 

.2072 

.3230 

■BOB 

.2634 

.3591 

.5500 

.0976 

.2392 

.3572 

mm  8 

.2964 

.4056 

.6000 

.0990 

.2436 

.3880 

.6500 

■(EES 

.3366 

.4517 

.6500 

.0974 

.2392 

.4213 

.  70?0 

.  1663 

.3753 

.4938 

.7000 

.0982 

.237? 

.4611 

.7500 

.  1807 

.4063 

.5308 

.7500 

.2430 

.5074 

.8000 

.1881 

.4261 

.5620 

.8000 

.1123 

.2571 

.5564 

.8500 

.1877 

.4330 

.5865 

.8500 

.2764 

.6024 

.9000 

.1797 

.4269 

.6031 

.9000 

.1324 

.2959 

.6410 

.9500 

.  1649 

.4089 

.6096 

.9500 

.1369 

.3093 

.6690 

1 .0000 

.1443 

.3810 

.6028 

1.0000 

.3109 

.6835 

1 .0500 

.1194 

.3453 

.5779 

1.0500 

.1207 

.2965 

.6802 

1 .1000 

.0915 

.3040 

.5270 

1.1000 

.0986 

.2647 

.6521 

1.1500 

.0616 

.2588 

.4355 

1 . 1500 

.0702 

.2184 

.5888 

1 .2000 

.0301 

.2105 

.2731 

1  .2000 

.0421 

.1659 

.484(1 

SIGMA* 

=  0.500 

SIGMA*  =  0.340 

.6000 

.0915 

.2226 

.3925 

.7000 

.0921 

.2345 

.4378 

.6500 

.  1016 

.2383 

.4411 

.7500 

.0970 

.2376 

.4872 

.7000 

.0995 

.2362 

.4691 

.8000 

.0919 

.2299 

.5089 

.7500 

.0945 

.2269 

.4972 

.8500 

.0861 

.2166 

.5378 

.8000 

.0920 

.2178 

.5379 

.9000 

.0839 

.2038 

.5899 

.8500 

.0938 

.2133 

.5951 

.9500 

.0857 

.  1965 

.6593 

.9000 

.0992 

.2149 

.6611 

1  .0000 

.0891 

.1972 

.7192 

.9500 

.1056 

.2215 

.7207 

1.0500 

.0902 

.2035 

.7404 

1 .0000 

.1096 

.2298 

.7591 

1  .  1000 

.0843 

.2067 

.7139 

1.0500 

.1078 

.2344 

.7683 

1.1500 

.0675 

.1899 

.6511 

1.1000 

.0974 

.2283 

.7466 

1.2000 

.0376 

.1258 

.5705 

1.1500 

.0773 

.2028 

.6970 

1.2000 

.0481 

.1483 

.6283 

*Cavitaf.ion  number. 
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TABLE  9  -  PROPELLER  OPERATING  POINTS 


Performance 

Parameters 

Propeller  Number 

Design 

4717B 

4717C 

4738A 

V,  knots 

59,0* 

59.0 

59.0 

59.0 

J 

1.037 

0.942 

0.982 

1.010 

n.  rpm 

1000 

1103 

1059 

1029 

n 

0.67 

0.52 

0.62 

0.72 

^Rounded  from  58. 

00 
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APPENDIX  A 

CAVITY  BOUNDARY  CONDITION 


From  the  Bernoulli  equation  with  respect  to  coordinates  fixed  on  the  blade 


(48) 


where  P  and  q  are  the  pressure  and  the  speed  of  a  fluid  particle  relative  to  the 
blade.  Thus 


<vO 


(<4+q- 


=  G  V 


A 


(49) 


Therefore 


AP 


2 G  \ 

V  V 
s  s 


(50) 


Again  from  the  Bernoulli  equation 


or 


P-P.  ,  o  2i 

—  -f  -it(V+uT)  +un) 


P 


2V  ^T 
V  V 


s  s 


(51) 


This  is  Equation  (25)  or  ('5) 


s;> 
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Writing 


Equation  (28)  can  be  written 


V£  = 


V  +  u, 


(m) 


+  u. 


(52) 


V 


s 


V 

s 


or 


(53) 


m 


is  as  small  as  or 


terms  including  G  u!£  and  G  uf 


can  be  neglected. 
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APPENDIX  B 
LEADING-EDGE  SOURCE 


We  assume  the  normal  velocity  on  the  foil  cavity  near  the  leading  edge  has  the 


a,  —  +  a„  +  a»x  =  b . 

1  /j  2  3  1 


solve  for  a^,  and  a^  with  three  points  of  x 


x  =  0.025,  0.05,  and  0.1 


i.  — — —  +  a,  +  a  0.025  =  b 

1  /OB' 


a,  - +  a„  +  a„  0.05  =  b„ 

1  /Os  2  3  2 


a,  •  -  •  ■  +  a„  +  a„  0.1  =  b_ 

1  /oa  2  3  3 


Then  with  the  determinant 


1  0.025 


'0. 025 


D  =  — ±2  1  0.05  =  0.0598745 

»05 


1  0.1 


we  obtain 
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a1  =  (0.05b1-0.075b2+p.025b3)/D 

a2  =  (-0.28911^+0. 5533985b2-0.2044243b3)/D 

a3  =  (1. 30985828^3. 1622776b2+1.8524193b3)/D 


The  total  flux  near  the  cavity  is 


o 


(V  -V  n )  dx  =  mx 
nu  nS/ 


rax  is  the  flux  of  the  corresponding  wedge  (or  uniformly  distributed  up  to  x) 
That  is 


m  = 


^2a^/x  +a2x+ 


2 


+  a2x+a 


This  is  the  same  as  the  leading  edge  cavity  thickness  at  x 
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